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Abstract 

A constant-amplitude code is a code that reduces the peak-to-average power ratio (PAPR) in 
multicode code-division multiple access (MC-CDMA) systems to the favorable value 1. In this paper 
quaternary constant-amplitude codes (codes over Z4) of length 2™ with error-correction capabilities are 
studied. These codes exist for every positive integer m, while binary constant-amplitude codes cannot 
exist if m is odd. Every word of such a code corresponds to a function from the binary m-tuples to Z4 
having the bent property, i.e., its Fourier transform has magnitudes 2™/^. Several constructions of such 
functions are presented, which are exploited in connection with algebraic codes over Z4 (in particular 
quaternary Reed-Muller, Kerdock, and Delsarte-Goethals codes) to construct families of quaternary 
constant-amplitude codes. Mappings from binary to quaternary constant-amplitude codes are presented 
as well. 

Index Terms 

Bent function, code, code-division multiple access (CDMA), Delsarte-Goethals, Kerdock, multi- 
code, peak-to-average power ratio (PAPR), quaternary, Reed-Muller 

I. Introduction 

Multicode code-division multiple access (MC-CDMA) is a simple scheme to implement rate 
adaption in CDMA systems [8]. The basic idea is to assign additional spreading sequences to a 
user who wishes to transmit with a higher data rate. In order to avoid self-interference, the user 
commonly employs n = 2™ orthogonal spreading sequences, which can be viewed as the rows of 
a Hadamard matrix of order n. Thus the data rate of this user is n times that in a conventional 
CDMA system. Typical values of m are 2 to 6. The principal drawback of this technique is 
that the transmitted signals can have a high peak-to-average power ratio (PAPR). That is, the 
peak transmit power can be much larger than the average transmit power. Thus the efficiency 
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of analog devices in the transmission chain, such as the power amplifier, digital-to-analog, and 
analog-to-digital converters, is limited due to the high PAPR of the signals. 

An elegant solution to solve this power-control problem is to draw the modulating words from 
a block code that contains only words with low PAPR and, simultaneously, has error-correction 
capabilities. This idea was originally proposed for orthogonal frequency-division multiplexing 
(OFDM) systems [9], where a similar power-control problem occurs. We will see that the PAPR 
is always at least 1. A code for which all codewords achieve this lower bound is called a 
constant-amplitude code. 

A coding-theoretic framework for binary codes in MC-CDMA has been established by Pa- 
terson [17]. It was shown that codewords with low PAPR are exactly those words that are 
far from the first-order Reed-Muller code, RM(l,m). This fact was used in [17] to prove 
fundamental bounds on the trade-off between PAPR, minimum distance, and rate of binary 
codes. Moreover several families of binary constant-amplitude codes have been constructed in 
[17] by exploiting the relation between bent functions [19], [13] and binary constant-amplitude 
codewords (a connection that was first recognized by Wada [24]). These codes are unions of 
cosets of RM(l,m) lying in higher-order Reed-Muller, Kerdock, or Delsarte-Goethals codes. 
Therefore they enjoy high minimum distance and are amenable to efficient encoding and decoding 
algorithms. 

Sole and Zinoviev [23] constructed binary codes with PAPR much greater than 1, which 
makes them less attractive for practical values of n. However in many situations their parameters 
asymptotically beat the Gilbert- Varshamov-style lower bound derived in [17]. 

While previous work was focused on binary codes, several motivations exist to study quater- 
nary codes for MC-CDMA. First, quaternary modulation rather than binary modulation is often 
employed in MC-CDMA systems [4]. Second, binary constant-amplitude codes cannot exist for 
lengths 2"* when m is odd (we will see that quaternary codes indeed do). 

In [20] the author established a connection between the coding problems in OFDM and in MC- 
CDMA. As a consequence, the generally nonbinary codes in [5], [16], [21], and [20] developed 
for OFDM can be re-used directly in MC-CDMA. However the best known upper bound on the 
PAPR of such codes is 2, and it can be shown using the proof of [20, Thm. 18] and the remarks 
thereafter that this approach cannot be used to construct quaternary constant-amplitude codes. 

The key concept in this paper is the connection between words with PAPR equal to 1 and 
functions from the binary m-tuples to Z4 having the bent property, i.e., the absolute values of 
their Fourier transform take on a constant value. This connection together with results from 
the theory of algebraic codes over Z4 is employed to construct families of quaternary constant- 
amplitude codes. Many of the resulting code families may be viewed as quaternary analogs of 
those developed in [17]. 
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The remainder of this paper is organized as follows. In Section HI] we will state the coding 
problem in MC-CDMA formally. In Section HlI] we establish our main notation and prove some 
basic properties of quaternary constant- amplitude codes. In Section ITVl mappings from binary to 
quaternary constant- amplitude codes are studied. In Sections |V] and |VI] we construct constant- 
amplitude codes from quaternary Reed-MuUer codes and from quaternary Kerdock and Delsarte- 
Goethals codes, respectively. Section IVIII contains some final remarks. 

II. System Model and Problem Statement 

We work with a simplified discrete-time MC-CDMA model, which essentially follows that in 
[17]. Let n := 2™, and let the canonical Walsh-Hadamard matrix of size n x n he recursively 
defined by 

Hn:=(y with H,:=il). (1) 

Let u; be a primitive 2''th root of unity in C, e.g., to = exp(z27r/2''), where = —1. Given a 
Z2h-valued word c = (cq, ci, . . . , c„_i), the transmitted MC-CDMA signal can be modeled as 

n-l 

S,{t) = J2^"'iHn)j,t, t = 0,1,..., n-l. (2) 

j=0 

In words, cj is used to modulate the jth row of Hn, and the transmitted signal is the sum of 
these modulated sequences. The PAPR of the word c (or of the signal Sc{t)) is defined as 

PAPR(c) := - max \SJt)\^. (3) 

n 0<t<n 

Observe that the PAPR is at least 1 and can be as much as n. 

A code C of length n over Z2/1 is defined as a subset C C Z^^. C is called linear if it is a 
subgroup of 7/^^ . We are particularly interested in the cases where /i = 1 or /i = 2, in which cases 
we say that C is binary or quaternary, respectively. If C is a binary code, then d//(C) denotes 
its minimum Hamming distance. Likewise, if C is a quaternary code, then dL(C) denotes its 
minimum Lee distance (see, e.g., [6] for details on Lee distance). The rate of C is defined as 

R{C) :=-log2|C|, 

n 

and the PAPR of C is defined to be 

PAPR(C) := maxPAPR(c). 

Our goal is to design quaternary codes C with PAPR(C) being much lower than n. Ideally 
PAPR(C) is equal to 1, in which case C is called constant-amplitude code. 
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III. Generalized Bent Functions 

Let 

F := {ze ■.z^ = z) 

be the set of TeichmuUer representatives in . Then every z E can be written uniquely in 
2-adic expansion 

h-l 

z = ^ Zj2^, where zo,zi,. . . , Zh^i e F. (4) 
i=o 

We define an operation on F by a © 6 := (a + 6)^'' \ Then (F, ©, ■) = (Z2, +, ■) is the binary 
field. The operation '+' is reserved to denote addition in Z2h. 

A generalized Boolean function is defined as a mapping / : F™ Z2h. Writing k = 
(fco, fci, • • . , fcm-i) for k G {0, 1}"*, every such function can be uniquely expressed in the 
polynomial form 

m— 1 

f{x) = f{xo,...,x^_i) = CfcJJx^', CfcGZah, 

fcG{0,l}™ j=0 

called the algebraic normal form of /. The degree of / is defined to be 

max vftnik), 

where wt^ik) is the Hamming weight of k. Functions having algebraic normal form 

where fco, fci, . . . , km-i G {0, 1}, are called monomials. Every generalized Boolean function can 
be expressed as a Z2h -linear combination of monomials. By ([H), every / : Z2h can be 

written uniquely in 2-adic expansion, viz 

h-i 

f{xo, . . . = ^/j(xo, . . . ,Xm_i)2^, (5) 

j=0 

where each fj is a mapping from F"^ to F. 

With each generalized Boolean function / : F"^ — > Z2h we associate a word of length 2™ 
with elements in Z2h. This word is obtained from the algebraic normal form of / by listing all 
the values f{x) as x ranges over F™ in lexicographic order. That is, if £ = Yl^=o ^j^-', where 
£0, ^1, • • • , £m-i e F, the Ah element of this word reads f{io, ii, . . . , im-i)- By convention, we 
shall denote the function and the associated word by the same symbol, and since there is a 
one-to-one correspondence between them, we often use the terms function and associated word 
interchangeably. 
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The rth-order Reed-Muller code RM(r, m) of length 2™ is defined as the binary linear code 
that is generated by the m-variate monomials of degree at most r [13, Ch. 13]. It contains 

codewords and has minimum Hamming distance 2™~^. 

The Fourier transform of a function / : — > is given by / : F'^ — > C with 



u-x 



where ■ denotes the scalar product in and is a primitive 2''th root of unity in C. The 
multiset {f{u) : u G F"^} is called the Fourier spectrum of /. 

Definition 3.1: A function f : F"^ ^ is bent if = 2™/2 for every u G F"^. 

The name bent function was coined by Rothaus [19]. His definition of bent functions applies 
to functions from F"^ to Z2. Kumar et al. [10] generalized this definition to functions from 
to Zg for arbitrary q. Another generalization of bent functions appeared in [14], where the 
functions are defined from F^ to {0, 1/2, 1, 3/2} = Z4 (as an additive group) and are called 
real-valued bent functions . The latter essentially coincide with Z4-valued bent functions defined 
above. 

The following theorem will be the key to obtain our code constructions later in this paper. 
Theorem 3.2: Let c : Z2fe be a generalized Boolean function. Then we have 

PAPR(c) = — max \c(u)\^. 

In particular the PAPR of c is equal to 1 if and only if c is a bent function. 

Proof: Observe that the elements of the Walsh-Hadamard matrix 7^2™, shown in ([T]), are 
given by 



where I = Yl^=o^ ^j'^^ ^^'^ ^ = Sj=o^ ^i^'' ^^^^ ^i' ^ ^ ^'^^ every j = 0,1, . . . ,m — 1. Thus, 
for each generalized Boolean function c : F™ Z2h, the corresponding MC-CDMA signal, 
given in satisfies 

Scit) = c{to, ti, . . . , trrt-l)- 

The first statement in the theorem then follows from ([3]). We immediately conclude that, PAPR(c) 
1 if c is a bent function. The converse is a consequence of Parseval's identity 

1 



2m / ■/ 



\c(u)\^ 



□ 
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We remark that in the binary case (where h = 1), Theorem 13.21 was essentially stated by Wada 
[24] and by Paterson [17, Thm. 8]. 

It is known (see, e.g., [19]) that a Z2-valued bent function cannot exist if m is odd and that 
the degree of a Z2-valued bent function is at most m/2 for m > 2. Therefore binary constant- 
amplitude codes can only exist for even m and must be contained in RM(m/2,m) if m > 2. 
This gives us an upper bound on the size of a binary constant amplitude code and, perhaps 
surprisingly, a lower bound on its minimum Hamming distance. In the remainder of this section 
we shall derive analogous bounds for quaternary constant- amplitude codes. 

Lemma 3.3: The values in the Fourier spectrum of a Z4-valued bent function on are of 
the form 2'^/^uj"'i^, where uj = {1 + i)/V2 and k G Z4. 

Proof: Let S denote an arbitrary value in the Fourier spectrum of a Z4-valued bent function 
on F"^. Then Re(5') and Im(5') must be integers and = 2™ must be a sum of two squares. 
From Jacobi's two-square theorem we know that 2*" has a unique representation as a sum of 
two squares, namely 2*" = (2'"/2)2 + 0^ if m is even, and 2™ = (2{"-i)/2)2 + (^2^m-i)/2y ^ 
is odd. Hence, if m is even, either Re(S') or Im(S') must be zero. If m is odd, we must have 
|Re(S')| = |Im(S')|, which proves the lemma. □ 

Recall from ([5]) that every / : F"^ Z4 can be written uniquely in 2-adic expansion, viz 

f{x) = a{x) + 2b{x), 

where a, b : F"^ — > F. We can express the Fourier transform of / in terms of the Fourier 
transforms of a and h as follows 



/(«)= i<-^+^'^-\-r 



u-x 



b{u) + (a © 6) (u) . b{u) - (a © 6) (u) 

H I ■ , (6) 



2 

which follows directly from the identity 

l + + (-1)1 



Theorem 3.4: With the notation as above, let m > 2 and suppose that / : F™ ^ Z4 is a bent 
function. Then the degrees of a and b are at most [m/2]. 

Proof: First suppose that m is even. Then, by Lemma 13. 3[ either the real part or the 
imaginary part of f{u) must be zero. This can only happen if b and a © 6 are bent, which 
implies that the degrees of a and b are at most m/2. Now let m be odd. By Lemma 13.31 in 
this case the real part and the imaginary part of f{u) must be equal in magnitude. This is only 
possible if the Fourier transforms of a © 6 and b take on only the values and ±2'^™"'"^)/^. Such 
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functions are called almost bent functions in [3]. It is known [3, Thm. 1] that the degree of such 
a function is at most (m + l)/2. □ 
Theorem 3.5: For m > 2 every quaternary constant-amplitude code Q of length 2™ satisfies 

r™/2i 



j=0 

and 

di(Q) >2L-/2J. 

Proof: By Theorem 13.41 every quaternary constant-amplitude code is contained inside the 

code 

{a + 26 : a,6 G RM([m/2],m)}, (7) 

whose minimum Lee distance is at least the minimum Hamming distance of RM([m/2],m). 
Since at least one word in the code ([7]) has PAPR greater than one (e.g., the all-zero codeword), 
it cannot be itself a constant-amplitude code. The theorem follows from the properties of 
RM([m/2],m). □ 
The preceding theorem provides a reasonable rate bound for small m. For instance, if m = 4, 
the bound asserts that the rate of a quaternary constant- amplitude is at most 21/16 (assuming 
an integer number of encoded bits). We will later construct a code with rate 18/16 (see Table H])- 
We speculate that the rate bound tends to be loose for large m. However, even for large m, the 
bound says that the constant-amplitude property of a quaternary code must be paid by a rate 
reduction from 2 to at most slightly more than 1. 

IV. Mappings From Binary to Quaternary Codes 

In this section we study mappings from binary to quaternary constant-amplitude codes, which 
can be used in connection with the results obtained in [17, Sec. V] to construct quaternary 
constant-amplitude codes. Recall from Theorem 13.21 that every constant-amplitude code must be 
comprised of bent functions. Therefore we shall study mappings from Z2-valued bent functions 
to Z4 -valued bent functions. 

Some mappings are based on the Gray map : Z4 ^ given by 

0^(00), 1^(01), 2^(11), 3^(10). 

In other words, writing z = a + 2h where a,b E F, we have 0(2;) = (6, a © 6). We extend 
naturally to act on words in Z4 . A useful property of (p is that it is a distance-preserving bijection 
from Z4 equipped with the Lee distance to F^" equipped with the Hamming distance [6]. That 
is, if Q is a quaternary code, we have 
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If / is a Z4-valued word of length 2™, we will often be interested in the Z2-valued Boolean 
function that describes </>(/). If we write f{x) in 2-adic expansion, viz /(x) = a(x) + 2b(x) 
where a, b are F-valued, then by the definition of the Gray map, 0(/) is a function from 
to F given by (j){f){x,y) = a{x)y © h{x). 

The theorem below has been essentially stated in [15, Thm. 2]. 

Theorem 4.1: Let m > 1, and let a, b : — ^ F be bent functions. Then / : Z4 
given by 

fix, y) = 2a(x)(l + y) + 2b{x)y + y 

is also bent. 

Proof: Similar to the proof of [15, Thm. 2]. □ 
Corollary 4.2: Let S be a binary constant-amplitude code of length 2™"^ (so m is necessarily 
odd). Then the code 

Q = {{2a + e ■ 1 , 2b + (1 + e) ■ 1) : a,b e B,e e F} 

of length 2™ is a quaternary constant-amplitude code with 

dL(Q) = 2d^^(i3) 
R{Q) = R{B) + 1/2™. 

Here, 1 denotes the all-one word of length 2™^^. 

Note that the quaternary code in the preceding corollary is essentially an offset binary code 
that is obtained by concatenating two suitably offset codewords of a binary constant-amplitude 
code. This has to be compared with [17, Lem. 19 and the remarks thereafter], where essentially 
the same code construction is suggested resulting in, however, codes with PAPR equal to 2. 

Theorem 4.3: Let a, b : — > F be bent functions, and define g : F™"*"^ — > F by 

g{x, y) = a{x)y © b{x){l © y). 

Then / = (f)~^{g) is a bent function. 
Proof: By rewriting g{x, y) as 

g{x, y) = b{x) © [a{x) © b{x)]y, 

it is seen that f{x) has 2-adic expansion 

/(x) = 2b{x) + [a{x) © 

Using we obtain 

, bin) + a(u) , - aiu) 
fw = + t . 
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Since a{u) and b{u) take on only the values ±2™/^, it follows that \f{u)\ = 2™/^ for each 

u e F™. □ 

The preceding theorem together with the distance-preserving property of the Gray map implies 
the following. 

Corollary 4.4: Let i3 be a binary constant-amplitude code of length 2™ (so m is necessarily 
even). Then the quaternary code 

Q = <l)-'{{{a,b):a,beB}) 

of length 2™ is a quaternary constant- amplitude code with 

dL{Q) = dniB) 
R{Q) = 2R{B). 

Theorem 4.5: Let m > 1, and let g : _p™+i — ;> F be a bent function. Then / = (p~^{g) is also 
bent. 

Proof: Recall that, when g is expressed as g{x, y) = a{x) y © b{x), where a, b : F"^ F, 
then / is given by /(x) = a(x) + 2b{x). The Fourier transform of g can be written as 

^(■U v) = ^_'^y(x)y®b{x)<Su-x<Svy 
xeF"',yeF 



^ ^ ^_'^-^b{x)(Bu-x _|_ j^_-|^^i;®a(x)-j 



Using the expansion ^ we obtain 

c- ?(m,0) , 1) 

Since, by assumption, \g{u,v) \ = 2^'"+^)/^ for each u E and each t> G F, we have \ f{u)\ = 
2^/2 fQj. g^^j^ ueF"^. □ 

With the distance-preserving property of the Gray map we obtain our next corollary. 
Corollary 4.6: Let i3 be a binary constant-amplitude code of length 2™+^ (so m must be odd). 
Then Q = (p~^{B) is a quaternary constant-amplitude code of length 2™ with 

dL{Q) = dniB) 
R{Q) = 2R{B). 



V. Codes From Quaternary Reed-Muller Codes 
Below we recall two quaternary generalizations of the binary Reed-Muller code from [5]. 
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Definition 5.1: 

a) For < r < m the code RM4(r, m) is defined as a linear code over Z4 of length 2"^ that is 
generated by the m-variate monomials of degree at most r [5]. 

b) For < r < m + 1 the code ZRM(r, m) is defined as a linear code over Z4 of length 2™ that 
is generated by the m-variate monomials of degree at most r — 1 together with two times 
the m-variate monomials of degree r with the convention that the monomials of degree — 1 
and m + 1 are equal to zero [5] (see also [6] for a slightly different, generally nonequivalent, 
definition of ZRM(r, m)). 

The code RM4(r, m) contains 

4^j=o j ) 

codewords and has minimum Lee distance 2"*"'', while ZRM(r, m) is a subcode of RM4(r, m) 
that contains 

codewords and has minimum Lee distance 2"*"^+^ [5]. 
The next lemma will be useful in the sequel. 

Lemma 5.2: Every word in a coset of ZRM(1, m) has the same PAPR. 

Proof: The lemma is a direct consequence of the fact that the absolute values in the Fourier 
spectrum of / : F"^ — > Z4 and the absolute values in the Fourier spectrum of a function given 
by 

f{x) + 2c-x + e, c e F"*, e e Z4 

are identical. □ 
The following theorem identifies a large number of Z4-valued bent functions. 
Theorem 5.3: Let cr be a permutation on F'', let gf : — > Z4 be arbitrary, and let / : F^^ — > 

Z4 be given by 

f{x,y) = 2a{x) - y + gi^x). 

Then / is a bent function. 
Proof: Write 

— i9{3:)^_i^u-x _ ^_'^\^{(T(x)+vyy 

xeF^ yeFl= 



X 



eF^ yeF^ 



11 



If X ^ V, the exponent in the inner sum is a linear form in y. Hence the inner sum is zero if 
X ^ V. On the other hand, if x = v, the sum is trivial and becomes equal to 2^. Thus 

and the proof is completed. □ 

The preceding theorem could be understood as a generalization of the Maiorana-McFarland 
construction of Z2-valued bent functions (see, e.g., [13, Ch. 14, Problem (20)]). A special case 
of this construction has been stated in [14, Thm. 1] in terms of real- valued bent functions. 

Now write the permutation a in Theorem 15.31 as ((Tq, 0"i7 • • • ? cTfc-i), where ao, . . . , au-i '■ — > 
F. Every aj must be balanced, so has degree at most k — 1 (see, e.g., [13, Ch. 13, Thm. 1]). 
Therefore the functions / identified in Theorem 15.31 have degree at most k. We deduce the 
following coding option. 

Construction 5.4: For m = 2k Theorem 15.31 identifies 2"^/^! ■ 42™'^^ bent functions, which 
organize in cosets of ZRM(l,m) inside RMi{m/2,m). The union of them is a constant- 
amplitude code of length 2™ with minimum Lee distance 2™/^. Assuming an integer number of 
encoded bits, for m = 4,6,8 the rate of this code is equal to 12/16, 31/64, and 76/256, 
respectively. By restricting the functions g in Theorem 15.31 such that the coefficient of the 
monomial of degree m/2 in the algebraic normal form of g is even, we can construct a subcode 
of the code considered above that is contained in ZRM(m/2,m). This code has half as many 
codewords as the code contained in RM4(m/2, m), but minimum Lee distance 2"*/^+^. The rate 
of this code is equal to 11/16, 30/64, and 75/256 for m = 4, 6, 8, respectively. 

Now we turn our attention to the code ZRM(2,m). This code comprises 2'^('"+^)/2 cosets of 
ZRM(1, m), where the coset representatives are given by a function Q : F^ Z4 that can be 
written as 

Q(x) = xBx'^ (8) 

m—l 

= ^^'J'^i + ^ ^jkXjXk, 

j=0 0<j<k<m 

where B = (bj^) is an F-valued symmetric matrix of size mxm. Such a Q is called a "Z^-valued 
quadratic form [1]. We say that Q has rank r and write rk(Q) = r if the matrix B has rank r 
(being computed over F). 

Theorem 5.5: With the notation as above, 

\Q{u)\ < 2™-^'^(«)/2_ 
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Proof: We have 



\ xgF'" ^ ^ y&F 

_ ^ j^xBx'^ -yBy'^ +2u-(x+y) 

x,y&F^ 

^ ^ ^(x+j/)B(x-j/)T'+2«-(x+y) 

Next we apply the variable substitution 

z := X ® y = X + y + 2x * y, 

where x * y denotes the component- wise product of x and y. Notice that for each fixed y, the 
mapping a; t-^ z is a bijection on F™. Using the identities 

x + y = z — 2x*y 

X — y = z — 2x*y — 2y, 

we can write 

= ^ ^iz+2x*y)B{z+2x*y+2y)T +2U-Z 
z,yeF"^ 

_ ^zBz'^+2zBy'^+2u-z 
z,y&F^ 

_ ^zBz'^+2u-z (^_\yBy'^ 

z&F"^ y£F"^ 

If the rank of B is r, then there exist 2'"^'' elements z E F"^ such that zBy"^ = (mod 2) 
for each y G and the inner sum is equal to 2™. On the other hand, for the remaining 2^ 
elements z G F™, the expression zBy-'^ is a nonzero linear form in y and the inner sum is zero. 
We conclude that 

|Q(u)p < 2"^ X 2"^-^ 

which proves the theorem. □ 
The preceding theorem in connection with Theorem 13. 21 and Lemma 15^ yields an upper bound 

on the PAPR of the cosets of ZRM(1, m) inside ZRM(2, m). 

Corollary 5.6: Let Q : ^ Z4 be a Z4-valued quadratic form. Then the PAPR of the coset 

Q + ZRM(l,m) is at most 2"'-'^'-Q\ 

A simple constant- amplitude code can now be constructed as follows. 

Construction 5. 7: Let Q : F™ — Z4 be a Z4-valued quadratic form of rank m. For instance, 
we may take 

Q{x) = Xo + Xi-\ h Xm-l 
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N{m) 



(in which case the corresponding B is the identity matrix). Then the coset Q + ZRM(l,m) is 
a constant-amplitude code that contains 2™+^ codewords and has minimum Lee distance 2"^. 

In what follows we will construct a constant-amplitude subcode of ZRM(2, m) by collecting all 
cosets of ZRM(l,m) inside ZRM(2,m) whose coset representatives correspond to Z4-valued 
quadratic forms of rank m. The number of such forms is equal to the number of m x m 
nonsingular symmetric matrices over F. This number is well known (see, e.g., [12]) and given 
by 

'm/2 

Y[(2"'+^ - 2^^) if mis even 

(m-l)/2 

JJ (2"^_22j) if mis odd. 
. j=o 

Construction 5.8: ZRM(2,m) contains N{m) cosets of ZRM(l,m) with PAPR equal to 1. 
The union of 2 1-^°^2 N(m)\ ^^.j^ cosets is a constant-amplitude code of length 2"* with minimum 
Lee distance 2"*~^ and rate 

Llog2 N{m)\ +m + 2 

2m 

For m = 4, 5, 6, 7 the rate is equal to 14/16, 20/32, 27/64, 35/128, respectively. A crude estimate 
yields N{m) > 2l-"*^/^J. Thus the rate of this code is bounded below by 

LmV2j +m + 2 

2m 

VI. Codes From Quaternary Kerdock and Delsarte-Goethals Codes 

We will begin this section with a short review of symmetric bilinear forms, which we will 
use to construct the quaternary Kerdock and Delsarte-Goethals codes. These constructions are 
then utilized to derive two classes of quaternary constant-amplitude codes. 

We let E be an extension field of F, so that {E, ©, ■) is isomorphic to GF(2"*). A symmetric 
bilinear form on £^ is a mapping B : E x E ^ F that satisfies symmetry, i.e., 

B(x,y) ^B(y,x) 

and the bilinearity condition 

B{x, ay ® bz) = aB{x, y) ® hB{x, z), Va, h e F. 

Note that bilinearity in the first argument follows from symmetry. The radical rad(i3) of B 
contains all elements x E E such that B{x, y) — for all y e It is a consequence of the 
bilinearity condition that rad(B) is a subspace of E. The rank of B is defined as the codimension 
of the radical, i.e., 

rk(B) :— m — dimrad(B). 
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Now let {Ao, Ai, . . . , A^_i} be a fixed basis for E over F. By the bilinearity condition we 
have 

(m—l m—1 \ m—1 m—1 

Vkh 1=00 ^M^^j, ^k)yk 
j=0 k=0 J j=0 k=0 

for all XQ.yo, . . . , x.rn-i,ym-i £ F. Therefore, relative to the fixed basis for E over F, the 
symmetric bilinear form B is completely determined by the m x m symmetric matrix B whose 
elements at positions (j, k) are given by B{Xj, A^). Hence we can write 

B{x, y) = (xo, . . . , Xm-i)B{yo, . . . , ym-if , 

where x = 0^1^^ XjXj and y = 0™Jo^ '^jVj- Note that the rank of the matrix B is equal to 
rk(i3). 

For the remainder of this section we define the linearized polynomial 

t 

La{x) := 0(%3;^' ® af ' x'^"' ') ® QqX, 

where a — (oq, Oi, ■ ■ ■ , a*) £ and <t < m/2. We also define the mapping Ba : E x E ^ 
F by 

^a{x,y) := tr(yL„(x)), 

where 

m—1 

tr(a;) := a;'' 
j=o 

is the trace function from £■ to F. It is not hard to verify that Ba obeys the conditions of a 
symmetric bilinear form. We define M.{t,m) to be the set of symmetric matrices that correspond 
to the set of bilinear forms {Ba : a e E^~^^}. The crucial property of M.{t,m) is given by the 
following theorem. 

Theorem 6.1: All nonzero matrices in M.{t, m) have rank at least m — 2t. 
Proof: Since Bq corresponds to the all-zero matrix, we assume a 7^ 0. Note that Ba{x, y) —0 
for every y e E if and only if La(x) — 0. Now observe that x 1— > x'* is an automorphism on E 
and 

t 

La{x^') = 0(aja;2*^' © aj™"'^''"') © agx'^' 

has degree at most 2'*. Therefore, La{x) has at most 2^* roots in E. We conclude that the 
number of x E E such that Ba{x,y) ~ for every |/ e -B is at most 2'*. Hence we have 
dimrad(,Ba) < It, so the rank of Ba (and therefore the rank of the corresponding matrix) is at 
least m — 2t. □ 
Next we use the set M.{t,m) to construct the quaternary Delsarte-Goethals code. 



15 



Definition 6.2: Let Q denote the set of Z4- valued quadratic forms corresponding to the 
symmetric matrices in A1(t,m) (cf. ([8])). For < t < m/2 the quaternary Delsarte-Goethals 
code is defined to be 

Vg{t,m) := [j g + ZRM(l,m). 

QeQ 

Moreover let /C(m) := VQ(0,m) be the quaternary Kerdock code. 

Our definition of VQ(t,m) essentially coincides with the definition given in [6], which we 
have extended from odd m to arbitrary m. We note that a similar definition of the quaternary 
Kerdock code also appears in [2]. 

The code VQ(t,m) is a linear subcode of ZRM(2,m) (because Ai{t,m) is closed under 
addition) and a union of 2'"(*+i) cosets of ZRM(1, m), so we have 

\Vg{t,m)\ = 2'"(*+2)+2. 

It is known (see, e.g., [7]) that 

dL(Vg{t,m)) = 2™-2*+L™/2J. 

By Theorem 16.11 and Corollary 15. 6[ we have the following. 

Corollary 6.3: The code Vg{t,m) \ ZRM(l,m) has PAPR at most 4*. 

In the particular case where t = we obtain a constant-amplitude code. 

Construction 6.4: The code /C(m) consists of 2™ cosets of ZRM(l,m), one of them is 
ZRM(l,m) itself. The union of 2™^^ remaining cosets is a constant-amplitude code of length 
2^" with rate (2m + l)/2'" and minimum Lee distance 2"* - 2L'"/2J. 

In the remainder of this section we shall construct a constant-amplitude subcode of VQil, m). 
The key to obtain the construction is given by the theorem below. 

Theorem 6.5: The number of nonsingular matrices in A^(l, m) is at least 4"*"^. 
Proof: For t = 1 we have 

La[x) = aix ® a^x ® x 

We are interested in the number of cases where the equation La{x) = has only the trivial 
solution X = 0. Squaring both sides and dividing both sides by x yields 

ajx^ © alx © ai = 0. (9) 

If ai = 0, there are 2"* — ! choices for ao such that La{x) = has exactly one solution, so assume 
in the following that ai 7^ 0. There are two cases. If ao = 0, ^ is equivalent to x^ = 1/ai. 
This equation has no solution if and only if a^^ ^ {x^ : x E E}, which can happen in 

0~gcd(3,2--l)) 
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cases. The latter expression is zero if m is odd and is equal to |(2™ — 1) if m is even. If ^ 0, 
(l9l) becomes under variable substitution x —x 

ai 

o a? 
x'^ + X + ^ = 0. 

% 

From [11, Appendix] we know that for each aQ E E \ {0} there exist |(2'" — (— 1)™) different 
ai E E \ {0} such that this equation has no solution. 

Straightforward counting yields that for either m the set Ai{l,m) contains 

2^ _ ^ ^ (2- + l)(2--l) ^ 

nonsingular matrices, which completes the proof. □ 
The preceding theorem leads to the following coding option. 

Construction 6.6: By Theorem 16. 5[ VQ{1, m) contains 2^"^"^ cosets of ZRM(1, m) that corre- 
spond to Z4-valued quadratic forms of rank m. The union of these cosets is a constant-amplitude 
code of length 2™, rate 3m/2™, and minimum Lee distance 2™ — 2^+L"i/2j 

Note added in proof: Recently, the author determined the complete rank distribution of the set 
M.{t,m) [22], which can be used to show that the number of nonsingular matrices in M.{t,m) 
is at least 2^*+^)™"^. Therefore, for general t, the code VQ(t,m) contains 

2{t+l)m-2 

cosets 

of ZRM(l,m) that correspond to Z4-valued quadratic forms of rank m. The union of these 

cosets is a constant- amplitude code of length 2"^, rate (t + 2)m/2™, and minimum Lee distance 

2»Tt 2*+L"^/2j 



VIL Final Remarks 

We have studied quaternary codes for R^PR reduction in MC-CDMA for the first time. Several 
constructions of codes with R^PR equal to 1 have been established. These codes are unions of 
cosets of ZRM(1, m). It is not hard to show that the image under the Gray map of such a code 
is a union of cosets of RM(1, m + 1). Since encoders and decoders for the binary image may be 
used to encode and decode, respectively, the corresponding quaternary code, existing encoding 
and decoding algorithms developed for binary OFDM codes [5], [18] are directly applicable 
here. 

For m E {4, 5, 6} the parameters of several constant-amplitude codes constructed in this 
paper are shown in Table IB This table illustrates how code rate can be traded against minimum 
Lee distance. Note that we have omitted those codes that were outperformed by one of the 
codes shown in the table. In particular the codes flowing from Construction 15.41 are notoriously 
surpassed by other codes; they appear to be rather weak at least for small m. 

Finally we wish to compare code rates and minimum Euclidean distances of our quaternary 
codes with those of previously constructed binary codes. Recall that the minimum Euclidean 
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TABLE I 

Parameters of Quaternary Constant-Amplitude Codes 



m 


rate 


min. Lee dist. 


reference 


4 


6/16 


16 


Construction |5J1 single coset of ZRM(1,4) 




9/16 


12 


Construction 16.41 subcode of K-l'i) 




14/16 


8 


Constructionist subcode of ZRM(2,4) 




18/16 


4 


[17, Construction 12] + Corollary 153 


5 


7/32 


32 


Construction 15.71 single coset of ZRM(1, 5) 




11/32 


28 


Construction 16.41 subcode of A;^(5) 




15/32 


24 


Construction 16.61 subcode of 'DQ{1, 5) 




20/32 


16 


Construction [5]8l subcode of ZRM(2, 5) 




23/32 


8 


[17, Construction 13] + Corollary |M] 


6 


8/64 


64 


Construction 15.71 single coset of ZRM(1,6) 




13/64 


56 


Construction 16.41 subcode of IC{6) 




18/64 


48 


Construction 16.61 subcode of T)Q{1, 6) 




27/64 


32 


Constructionist subcode of ZRM(2, 6) 




30/64 


24 


[17, Construction 18] + Corollary |4H 




40/64 


16 


[17, Construction 12] + Corollary 1441 




46/64 


8 


[17, Construction 13] + Corollary I4H 



distance of a binary (quaternary) code is equal to four (two) times its minimum Hamming (Lee) 
distance. It can be observed that all binary constant-amplitude codes proposed in [17, Sec. V] 
are outperformed by a quaternary constant-amplitude code constructed in the present paper. For 
example, [17] reports a code of length 64 with minimum Euclidean distance 32 and rate 23/64, 
while Table U contains a quaternary code with the same length and minimum Euclidean distance 
but rate 40 /64. 
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